YITP-95-10 



in 

On 
OS 



O 



in 

Os 

cr 



X 



Back- Reaction on the Topological Degrees 
of Freedom in (2+l)-Dimensional Spacetime 



Masafumi Seriu 



O \ Inter-University Centre for Astronomy and Astropkys.es 

Post Bag 4, Ganeshkhind, Pune 411001, India 

and 

> ' 

^ ■ Yukawa Institute for Theoretical Physics 



Kyoto University, Kyoto 606, Japan 



To appear in Physical Review D. 



* Present address. 



ABSTRACT 



We investigate the back-reaction effect of the quantum field on the topological 
degrees of freedom in (2+l)-dimensional toroidal universe, Ai — T 2 x R. Con- 
structing a homogeneous model of the toroidal universe, we examine explicitly 
the back-reaction effect of the Casimir energy of a massless, conformally coupled 
scalar field, with a conformal vacuum. The back-reaction causes an instability of 
the universe: The torus becomes thinner and thinner as it evolves, while its total 
2-volume (area) becomes smaller and smaller. The back-reaction caused by the 
Casimir energy can be compared with the influence of the negative cosmological 
constant: Both of them make the system unstable and the torus becomes thinner 
and thinner in shape. On the other hand, the Casimir energy is a complicated func- 
tion of the Teichmuller parameters (r^r 2 ) causing highly non-trivial dynamical 
evolutions, while the cosmological constant is simply a constant. 

Since the spatial section is a 2-torus, we shall write down the partition func- 
tion of this system, fixing the path-integral measure for gravity modes, with the 
help of the techniques developed in string theories. We show explicitly that the 
partition function expressed in terms of the canonical variables corresponding to 
the (redundantly large) original phase space, is reduced to the partition function 
defined in terms of the physical-phase-space variables with a standard Liouville 
measure. This result is compatible with the general theory of the path integral for 
the lst-class constrained systems. 
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1. Introduction 



Topological considerations are necessary in many situations. Since physical 
laws are usually expressed in terms of local, differential equations, their impor- 
tance is not prominent at first sight. However, once one proceeds to solve them, 
one has to take boundary conditions into account, which allow the topological 
information to enter in the theory. In general relativity, which handles the dynam- 
ics of spacetime, the topological properties acquire dynamical meaning and their 
consideration becomes more significant. The aim of this paper is to present an 
explicit, detailed investigation of the dynamics of topological degrees of freedom in 
spacetime, in the context of the back-reaction problem in semiclassical gravity. We 
concentrate on the case of (2 + l)-dimensional toroidal spacetime, M. ~ T 2 x R, 
and make use of various techniques developed for the 2- and 3-dimensional grav- 
ity. Here, we do not discuss the topology change [1][14]. The term "topological 
degrees of freedom", indicates those global parameters, describing the global de- 
formations of the spatial hypersurface, which are of topological origin (the moduli 
deformations) (§§3 — b and Appendix A). 

As a first preliminary study for the full quantum gravity, it is reasonable to 
consider the effect of the curvature of a fixed background spacetime on the behavior 
of quantum matter field, which is the subject of quantum field theory on a curved 
spacetime [2] [3]. Then, the next natural step is the investigation of the influence 
of such a quantum field on classical spacetime geometry, which is called the back- 
reaction problem in semiclassical gravity. Usually, one tries to describe this effect 
by the semiclassical Einstein equation, 

Ga[3 = OL < T a p > , (1) 

where < T a p > is some c-number, obtained from the energy-momentum tensor 
operator and the inner-product of some quantum states, and a is an appropriate 
gravitational constant with physical dimension [a] = [(length)" -2 ]. (Here, n is 
the spacetime dimension. We treat H as [h] = [1], and set c = 1 in this paper.) 
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There are several uncertain issues and technically complicated points about this 
treatment. First, it is not clear what kind of quantity should be chosen for < T a p > 
[4]. Here, we regard that < T a p > should be some expectation value, rather 
than the quantity < out\T a p\in >, since the latter harms the reality and the 
causal nature of eq.(l) [5] [6] [7]. Then, if one regards the path-integral formalism as 
fundamental for quantum gravity, the so-called in-in formalism [8] [5] [6] should be of 
more importance than the standard in-out formalism [7] . Second, the regularization 
of < T a p > requires complicated, though well-established, techniques, which itself 
is one main topic of the quantum field theory on a curved spacetime [2] [3] . Third, 
eq.(l) in general becomes complicated, even though < T a p > has been successfully 
computed, so that it is difficult to solve it and study the effect of the back-reaction 
in detail. Fourth, one can show that eq.(l) can be obtained from the first variation 
of the phase part in the in-in path-integral expression [5] [7] , in which matter part 
has been integrated out formally while gravity part is left unintegrated without 
the explicit fixation of the measure. If one wants to go one step further, however, 
one should also take care of the effect coming from the path-integral measure for 
the gravity part. It is usually difficult since a reasonable, general measure has not 
been fixed yet. Fifth, to speak rigorously, eq.(l) itself contains an inconsistency 
from the very beginning. Since gravity and matter couple, quantum fluctuations 
of matter cause corresponding quantum fluctuations of gravity. Thus, there is a 
limitation in principle to the semiclassical treatment (eq.(l)), because we try to 
treat gravity classically while matter is treated by quantum theory [7]. Specifying 
the exact validity conditions for eq.(l) is one of the main topics of semiclassical 
gravity [7] [9] [10]. 

In this paper, we consider a (2 + l)-dimensional spacetime M. ~ X x R, with 
E ~ T 2 , a torus. We choose, as a matter field, a massless conformally coupled 
scalar field with a conformal vacuum, and investigate explicitly the back-reaction 
effect, resulting from the Casimir energy of matter, on the topological degrees of 
freedom, i.e. the modular-deformations of the torus. As is stated above, the topo- 
logical degrees of freedom is one of the essential ingredients of spacetime dynamics. 
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However, the back-reaction on topological modes has seldom been discussed so far, 
partially because such a finite number of degrees of freedom are hidden in infinite 
number of gravity modes in 4-dimensional spacetime. One advantage of the re- 
duction of the number of dimension from 4 to 3 is that, only the finite topological 
modes plus a spatial volume remain dynamical for the case of pure gravity, due to 
the dimensionality [11] [12] [13]. One can understand this point as follows: When 
n — 3, the spatial metric has 3 independent components at each spatial point, 
while there are 3 constraints at each point. Thus, redundant infinite number of 
modes are gauged away and only a finite number of modes remains. Here, we want 
to investigate the back-reaction effect from matter onto the topological degrees of 
freedom of spacetime, which would force us to take the matter field into account. 
To preserve the above-mentioned nice property of the finiteness of the number of 
degrees of freedom, we choose a model in which the matter field is in a vacuum state 
on a spatially homogeneous (2+l)-dimensional spacetime. Another advantage of 
the reduction of dimension in the discussion of topological aspects comes from the 
fact that 2-dimensional topology is completely classified in a simple manner so that 
it is easy to construct various topologies [14]. 

Another good point of this model is that some difficulties and complications 
stated above of the semiclassical Einstein equation, eq.(l), become simplified and 
tractable to a great extent in this case: 

First, we choose a conformal vacuum |0 > as a natural candidate for a vacuum 
state of matter in our case, and use < 0|T a/ 3|0 > on the right-hand side of eq.(l). 

Second, since the background spacetime shall be chosen as (conformally) flat 
and the matter field is conformally invariant, < T a p{g) > can be calculated from 

< T a p{rj) > (rj: a flat metric) along with the trace-anomaly [2], which which 
simplifies the manipulation. Furthermore, in our case, the spacetime dimension is 
odd, n = 3, so that there is no trace-anomaly [2]. Thus, < T a p{g) > is related to 

< T a p{r]) > in a simple manner. 

Third, because of the dimensionality, eq.(l) is reduced to a set of six first-order 
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ordinary differential equations and we can investigate the effect of the back-reaction 
explicitly. 

Fourth, we restrict the metrics to a special class, with spatial part being the 
one for the locally flat metrics on a torus. Thus, we can fix the path-integral 
measure by the use of the techniques developed in string theories [15] [16]. Within 
this model, we can discuss explicitly the influence on the semiclassical dynamics of 
gravity. Our treatment corresponds to the minisuperspace approach in quantum 
cosmology: Putting restrictions on the variables to be quantized (e.g. spatial ho- 
mogeneity), which is compatible with the classical dynamics, quantum theory is to 
be constructed within this restricted subclass of variables. Though this treatment 
is self-consistent as a quantum system, one significant question naturally arises: To 
what extent such a treatment reflects faithfully the original full quantum theory? 
From the viewpoint of the original full system, the restrictions are regarded as 
constraints on the phase space, which can modify the path integral measure for 
the reduced variables (minisuperspace variables). Our model may be a good test 
candidate to investigate this point in detail. 

Fifth, the (in-in) effective action for gravity, W[g + : <?_], becomes relatively 
simple in our case, and this reduces to W[t+, 7+, V+ : rl, t2, V_], a functional of 
six functions of t, (t±, t±, V±), where V± indicate the spatial 2-volume (area) and 
(t±,t±) are the Teichmiiller parameters describing the topological degrees of free- 
dom of a torus. Although the exact calculation of W has already become difficult, 
we can still estimate its functional form to leading order in Ti. In computing W, 
our model reveals explicitly the peculiarity of the semiclassical gravity, compared 
with the standard treatment of the quantum dissipative system, e.g. the Brown- 
ian motion [17]: There is no linear coupling between the sub-system (gravity) and 
the environment (matter field). Their coupling is put in the kinetic term of the 
matter field. This model might provide the simplest non-trivial example for the 
investigation of the quantum dissipative system including gravity. 

In §2, we recapitulate how to handle quantum fields on topologically non-trivial 
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spaces: Construct the quantum field theory on M. ~ T 2 x R and calculate the 
Casimir energy of a massless, conformally coupled scalar field with a conformal 
vacuum. [2], [3], [19] 

In §3, we extract explicitly the topological degrees of freedom of a torus and 
reduce eq.(l) to a canonical system with a finite number of degrees of freedom 
[13]. Then, we investigate explicitly the effect of the back-reaction of matter on 
the dynamics of the topological degrees of freedom. We shall see that the back- 
reaction makes the system unstable and the torus becomes thinner and thinner as 
it evolves, while its 2-volume becomes smaller and smaller. These behaviors are 
universal that is independent of the initial conditions. The asymptotic analysis of 
the set of dynamical equations justifies this point. We shall also compare our case 
of the Casimir energy with the case of the negative cosmological constant, since 
both of them can be regarded as negative energies. Most significant difference is 
that the Casimir energy is a complicated function of the Teichmuller parameters 
(t^t 2 ), while the negative cosmological constant is just a constant. 

In §4, we investigate the partition function of this system, fixing the mea- 
sure with the help of the techniques in string theories. We show explicitly that 
the gauge-fixing reduces the partition function formally expressed in terms of the 
canonical variables for the (redundantly large) original phase space, to the parti- 
tion function defined in terms of the physical-phase-space variables with a standard 
Liouville measure. This result is compatible with the general theory of the path 
integral for the lst-class constrained systems. We also estimate the functional form 
of W to leading order in H. Section 5 is reserved for discussions. 
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2. Quantum field theory on a (2+l)-dimensional toroidal spacetime 

This section is for defining the model to be considered, and calculating the 
energy-momentum tensor in our model, as a preliminary for the next section, where 
the back-reaction effect is analyzed in detail. Calculating < T a p > is now a well- 
established topic, and we just sketch the essence in the context of our model for 
later uses. 

(a) Scalar field on a torus 

We consider a (2+l)-dimensional spacetime with topology T 2 x R. We con- 
centrate on the case when the geometry of the space E ~ T 2 is locally flat. A flat 
2- geometry is endowed on E by giving a metric* , 

dl 2 = h ah d£ a d£ h , (2 - a) 

where 




(2-6) 



and the periodicities in the coordinates £ x and £ 2 with period 1 are understood. 

c 

Here , (r^r 2 ) are the Teichmiiller parameters [15] [16] independent of spatial 
coordinates (£\£ 2 ), and r := r 1 + ir 2 , r 2 > 0. Note that Vh := (&eth ah ) 1 / 2 = 1. 

The Laplacian operator A := — 1/ \fh.d a {h ab \/lidi ) ) on E with the line element 
dl 2 (eq.(2 — a, b)) gives the normalized eigenfunctions 

/mn 2 (0 = exp(i27rni^ 1 ) • exp(?27m 2 £ 2 ) (ni,n 2 e Z) , (4) 

* For dcfinitcness, wc shall use the symbol h a b to represent the particular matrix given by 
(2 — 6), while the symbol h a b shall be reserved for more general context, representing a 
general spatial metric induced on a spatial surface X. 

§ Throughout this paper, r 2 always indicates the second component of (t^t 2 ), and not the 
square of r. The latter never appears in the formulae. 
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and the eigenvalues 



47T 2 

A ni n 2 = -^(k| 2 n 2 - 2r 1 nin 2 + n 2 ) . (5) 



Now, let us consider a spacetime M. ~ £ x R, with a line element <is 2 = 
—dt 2 +h a ( ) dC, a dC, b . The fundamental positive frequency solutions for Du(t, £ 2 ) = 
are 

^(«,0 = ^e-^«/ A (0 , (6) 

where A stands for nin2 and cj^ := V%1 — \f^nin 2 - Afterwards, we follow the 
standard procedure for the field quantization [2] [3] . 

(b) The model 

We shall investigate the back-reaction of the matter field on the topological de- 
grees of freedom (r 1 , r 2 ). The most ideal treatment of the back-reaction described 
by eq.(l) may be the self-consistent determination of the geometry g a p through 
eq.(l): < T a p > depends on g a p, and this g a p is self-consistently determined by 
eq.(l). However, it turns out that such a treatment becomes highly complicated 
even in our simple model. To make our analysis tractable, then, we treat the 
back-reaction in the following sense, which is usually adopted in the back-reaction 
problems [2] [3]: We prepare a background spacetime and calculate < T a p > on it. 
Then, we discuss the modification of the background geometry due to the < T a p >, 
using eq.(l). 

Now, as a background spacetime, we choose a solution of the vacuum Einstein 
equation, G a p = 0. More specifically, we prepare a locally flat spacetime, ds 2 = 
-dt' 2 + Vdl 2 = V(-dt 2 + dl 2 ), where dl 2 is given by eqs.(2-a, b), and V, r 1 and r 2 
are chosen to be constant for the background spacetime. (Below, we occasionally 

* In connection with the later applications, it is worthwhile to note that, even though r 
would depend on t, the form of the equation Dip — would not change, because of the form 
of the metric, g a p = (— 1, h a b) with detg a p = — 1. 
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treat this flat spacetime as conformally flat, just for mathematical convenience.) 
We choose as a matter field, a massless conformally coupled scalar field ip, 

S m = -\J (g a/3 d a ijd p ij + ^Rij 2 )V=gd 3 x . (7) 

The (improved) energy-momentum tensor operator [3] becomes, 
T a p{g) = ^^c^-^d^d 7 ^^ 

(8) 

We choose the conformal vacuum as a vacuum state for the matter field. Then, 
< T a p(g) > is simply related to < T a p(r)) > as, 

< T aP {g) >= V- 1 ' 2 < T ap (rj) > , (9) 

when the metric g a p and the flat metric r) a p are related as g a p = Vr] a p^ On flat 
spacetime, the field equation for ip becomes O0 = 0, and eq.(6) can be used as 
fundamental solutions. In this manner, the time evolution of V causes no direct 
complication in the analysis. 

However, the time-dependence of (r 1 ^ 2 ) caused by the back-reaction harms 
the self-consistency of the analysis, which is inevitable if the tractability of the 
back-reaction problem, described by eq.(l), is to be maintained. When (r^r 2 ) 
evolve in time, the functions in eq.(6) are no longer exact solutions for Otp = 0, 
because uja '■= V%1 becomes t-dependent, through the t-dependence of (r^r 2 ) 
(eq.(7)). Furthermore, the spacetime described by ds 2 = —dt 2 + V(t)h a bd^ a d^ b is 
no longer conformally flat when (r 1 , r 2 ) evolves, because of the t-dependence of h a b. 
Thus, we should look at the results of the analysis in an adiabatic sense, i.e. valid 
when terms including f 1 and f 2 are not dominant in the formulae prominently. 

* This simplification occurs, because < T a p{g) > for a conformally invariant field, with the 
conformal vacuum, on a conformally flat spacetime is completely determined by < T a p(r)) > 
and the trace anomaly < T£(g) >, while the latter vanishes when the spacetime dimension 
is odd [2]. 
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Such a conflict between self-consistency and the tractability of the analysis always 
occurs in the back-reaction problem. In our present model, this adiabatic treatment 
provides a good approximation because f 1 and f 2 , caused by the back- reaction, 
turn out to be sufficiently small (see §§3 — c). 

We next need Hadamard's elementary function [2] [3] G^\x) for ds 2 = —dt 2 + 
dl 2 to calculate < T a p(rj) >. This function and the related energy- momentum 
tensor have already been extensively investigated [19]. We first compute G^\x) 
for A4 ~ R 3 , and afterwards take care of the periodicity in A4 ~T 2 xR, adding all 
contributions from points which should be identified [18] [19]. For the 3-dimensional 
Minkowski space, G (i) 

(x) is, 

GW(x):=<0\{^(x),^(y)}\0>=^-(2a)- 1 / 2 (a > 0) , (10) 
where a := \x 2 = \r\ a ^x a x^ ', \ times a square of a world distance. Thus we get* 



t 00 ' 

G(1) ^) = ^ E (2<w*)r 1/2 , (ii) 

711,712=— 00 



where 



2<7 nin2 (x) := -t 2 + -1 |(^ + m) + r(£ 2 + n 2 )f 



Now it is straightforward to compute < T a p(r)) > explicitly. {r] a p = (— l,/i a &) 
with (2-6).) The result is, 

<roo> ._^ E ^^ , (12 _ a) 

nin2 1 1 



* The prime attached to the S-symbol, like in cq.(ll), indicates that the zero-mode (ni = 
n 2 = 0) should be excluded from the summation, whenever it causes a divergence. 

§ For computation it is helpful to note that < T a p(ri) >= \d a df3>Gty, where d a dp>Gty := 
d x .d xie G^{x-x') w ^ x andx" := (i,£\£ 2 ). 
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<ru> ^^! Ei i -M E (;^<f , (12 _ 6) 

An ^-^ \n\ + TU2 r An \n\ + TU2 \° 



n\ri2 niW2 



(tV + 


T 


w 




ni + ™ 2 


5 



(ni + r 1 


n2)(r 1 m + 


T 


2 n 2 ) 




ni + rn2 
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<T22 . ^W y 1 S^r 2 ) 1 / 2 ^ 

47r ^-^ Ini+rnol 3 47r ^-^ 

niri2 nin 2 

(12 -c) 

<T 12 >=< T 21 > 

An ^ \n\ + TU2\ S An ^ 

niri2 ' 1 nin 2 

(12 - d) 

< T 0a >=< T aQ >= (a = 1, 2) . (12 - e) 

For a metric g a/ 3 = V(— 1, h a b), < T a p{g) > can be obtained by eq.(9). Since the 
Planck scale is the only scale which comes into our model, we understand that a 
suitable power of a : = /pi an ck is multiplied to quantities like those in eqs.(12 — 
a) — (12 — e), if necessary, in order to adjust their physical dimensions. These 
contributions of order h to < T a p > in eqs.(12 — a) — (12 — e) originate from a 
non-trivial spatial topology £ ~ T 2 , and are well-known as the Casimir effect [2] [3]. 



3. Back-reaction of the Casimir effect 
on the topological degrees of freedom 

(a) The extraction of dynamics of the modular deformations 

Having computed < T a p(g) > in the previous section, we shall now investigate 
the back-reaction of < T a p(g) > on the evolution of the spacetime. We consider 
the Einstein gravity on M. ~ T 2 x R and a massless conformally coupled scalar 
field on it; S — ^ J R^J—g + S m , where a := /pi an ck an d S m is given by eq. 
(7). The canonical formulation is suitable to investigate the temporal evolution 
of the spacetime. We thus perform the (2+l)-decomposition, but care should 
be taken because of the presence of the conformally coupled field. In the back- 
reaction problem, we regard that ip 2 (x) is replaced by a vacuum expectation value 
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< vp 2 (x) >, which is independent of spatial coordinates. Furthermore, we shall 
finally choose the spatial coordinates s.t. iV a = so that n a = (— 1/N, 0). These 
facts simplify the procedure of (2 + l)-decomposition. 

Following the standard manipulation [20], we finally get the total action in 
canonical form, 

S = J 7i ab h ab -NH- N a H a , (13 - a) 

where the Hamiltonian constraint and the momentum constraint become, respec- 
tively, 

H = {{K ah K ab - K 2 -i 2 )R)/ a + <T aP > n a vP}^fh , (13-6) 
n a /Vh = -2D b (K a b -5 a b K)/a-<T aP >n (3 . (13 - c) 

Here, N, N a are the lapse and the shift functions, n a = (— 1/N, N a /N) is the 
normal unit vector of the spatial surface, and stands for the scalar curvature 
for the spatial surface E. The operator D a is the covariant derivative w.r.t. h a b, 
and 7r ab := (K ab — Kh ab )y/h/a, K ab is the extrinsic curvature of a spatial surface* 

We choose a coordinate system s.t. iV a = so that n a = (— 1/N, 0). Thus, 

< T a p > rfi = —1/N. < T a Q >= (a = 1,2) from eq.(12 — e). In our case, thus, 
the momentum constraint becomes 

HjVh = -2D b (K a b - 5 a b K)/a = . (13 - c/) 

Then, we can extract the moduli degrees of freedom (corresponding to the global 
deformations of a torus) by solving eq.(13 — c') explicitly [13]. 

* Throughout this paper, we use a spatial metric h a b, an induced metric on a spatial surface 
E, to raise and lower the spatial indices, a, b, c, and to define the spatial covariant 
derivative D a . In particular, the geometry of our concern is given by the line element, 
ds 2 = —dt 2 + Vdl 2 , with (2 — a, b). Thus, the spatial metric in our model is h a b = Vh a b, 
with (2-6). 
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The system of coordinates in our model (ds 2 = —dt 2 + Vdl 2 with (2 — a, b)) 
corresponds to York's time-slicing [21], i.e. the time-slicing by the spatial surfaces 
on which 

a := —K/a = const . (14) 
Thus, eq.(13 — c') is equivalent to 

an a /Vh = -2D b K b = , (13 - c") 

where K b := K b -\b b K, the traceless part of K b . It means that* K ab G KerP^, 
so that K ab can be expanded in terms of the basis of KerP\\ {& Aab }A=i,2 ! 

K ab = -j2PA^ ab ■ (15) 

In our case, E ~ T 2 , we can choose the lapse function as iV = N(t) without 
any contradiction with the York's slice. This is shown almost in the same manner 
as for the case of pure (2 + l)-gravity [12] [13]. Now, using some basic facts on the 
moduli space M.g=\ (see Appendix A), it is straightforward [13] to show that our 
system is reduced to 



n f , dV ^ ^ 

2 

_ M( ST g AB p APB - \a 2 o 2 V + a<T a(5 > n a n?V). (16) 
A,B=\ 

Note that the contribution from the spatial diffeomorphism has been eliminated 
from dynamics, by solving the momentum constraint (13 — c') explicitly. Only the 
Weyl deformations and the modular deformations have remained. 



★ See Appendix A for the terminology and notations related to the moduli space. 
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(b) The evolution of the Teichmiiller parameters caused by the back- 
reaction 

In our model, ds 2 = V{-dt' 2 + h ab d£ a d£ b ). Thus, from eq.(9) and n a = 
(-1/VV,0), we get < T al3 > n a n^ = V^l 2 < T 00 >, where < T 00 > is given by 
eq.(12 — a). (Note that this combination is coordinate independent.) 

By setting N(t) = 1, we get the canonical equations of motion described by 
the constraint function, 

aH = £ g AB PAPB - \a 2 a 2 V - ha{r 2 fl 2 f{r)V^I 2 = , (17) 



A,B=1 



where 



oo, 



47r ^-^ Ini + rnol 3 

711,712=— OO 

oo, 

1 V 1 

~^ ni h-oo (n? + 2rinin2 + |r|2^)3/2 " 

Clearly, /(-r\r 2 ) = /(r\r 2 ), /(r^n.r 2 ) = /(t\t 2 ), /(n+a,r 2 ) = /(n-a,r 2 ) 
(n: integer, a: real) and f(r l ,T 2 ) is singular at (t^t 2 ) = (n, 0). Furthermore, 
the combination 27r(r 2 ) 3 / 2 /(r 1 , r 2 ) appearing in (17) is equivalent to the non- 
holomorphic Eisenstein series G(r, 3/2), whose modular invariance as well as other 
properties are well-known [22]. The first term in eq.(17) is also modular invariant, 
since it behaves as a scalar field on the moduli space* Thus, the Hamiltonian 
constraint eq.(17) is modular invariant as it should be. Figures 1 — a, b show the 
behavior of the function /(t^t 2 ). 



* Another convenient way for discussing the invariance is to perform the Legcndre trans- 
formation of the action in concern, and to look at the action in terms of the configuration 
variables. In this case, the kinetic term for (r 1 ^ 2 ) becomes proportional to ^2Gabt a t b , 
which is clearly modular invariant. For the discussion in the context of the path- integral, 
including the discussion on the path-integral measure, see §§4 — a. 
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For the explicit investigation of the dynamics, let us first calculate G AB ac- 

v ( 1 r '\ 

cording to eqs.(A5) and (A2 — c) with = I I. (Note that 



T T 



rfs 2 = -dt 2 + Vdl 2 .) Then, we get 



v fo 1 \ v ( -i -t 1 

T ^=^2^ 2r lJ' T ^=^p)2^_ r l (r2 ) 2 - (r l) 2 J • (19 " a) 

Note that {TAab}A=i,2 are symmetric, traceless 2-tensors satisfying —ID^Aa — 
—2dt,TAa = 0. Thus, {TAab\ A=\,2 can also be utilized to form a basis for KerP\\ 
{$> Aab } A=l,2- By normalizing them to satisfy (^ a ,Tb) = S A b, we obtain, 



i _ 1 ( r M 2 _ i ( - 1 - rl , 

*«*-2^ T 2 2rVj' (r 2 ) 2 - (r 1 ) 2 ' ' (19 " &) 

Thus, the Weil-Peterson metric reduces to the one which is conformally equivalent 
to the Poincare metric, 

2V fl 

Gab = (7k, 7b) = 9 9 9 



(20) 



Hence, the geometry conformal to the Poincare geometry [23] (negative constant 
curvature geometry) is endowed on the Teichmuller space, which is equivalent to 
the upper half-plane H+ ((r 1 , r 2 ) G R x R + ). Then, the system has been finally 
reduced to the constraint system ((V, a), (t 1 ,;?!), (r 2 ,p2); H = 0) with 

«H = ^^{v\ + V\) - \a 2 a 2 V - ha(r 2 f/ 2 f{r)V^ 2 = . (21) 
The equations of motion for (V, a) are 

V = -aaV (22 - a) 



1G 




(p? + ^)-f(r 2 ) 3/2 /(r)V- 3 / 2 



(22 - b) 



The equations of motion for (t 1 ,^) and (r 2 ,^) are 



a 




(23 - a) 



V 



df(r) 




(23 - b) 




(24 - a) 



P2 = —ApI + pI) + f (r 2 )V2 / (r)y-V2 + ^(,2)3/2^^-1/2(24 _ b) 



First, we should note that the time evolution becomes trivial when there is no 
matter field, /(r) = 0, in the following sense: In this case, eqs.(22 — a, b) allow 
a solution, a = 0, V — const, p\ = p2 = 0. It is clear that, from eqs.(21), 
(23 — a, b) and (24 — a, b), equations of motion do not allow any solution, compatible 
with cr = 0, V = const, other than r 1 =const, t 2 =const. This corresponds 
to the 3-dimensional Minkowski space in the standard coordinates (T, X^X 2 ) 
with suitable identifications in spatial section (X^X 2 ) described by (r^r 2 ). The 
unique solution above shows that there is no time evolution with respect to the 
standard time-slice, T =const ((7 = 0). This configuration is what we have chosen 
as a background spacetime. (However, there are different solutions characterized 
by the initial condition a ^ 0. In these cases, (r 1 , r 2 ) evolve in time.) 

The back-reaction of the quantum field causes a non-trivial evolution of (r 1 , r 2 ), 
i.e. global deformations of a torus. It is clear from eq.(21) that even when a m 
so that the term — ^a 2 <j 2 V in eq.(21) can be neglected, a non-trivial evolution 
of [t 1 ,t 2 ) occurs, because of the negativity of the term — /ia(r 2 ) 3 / 2 /(r)l /_1 ' /2 in 
eq.(21). The choice of the solution (7 = 0, V = const is not allowed any more, as 
is seen from eqs.(22 — a,b). 

Figures 2 — a, 6, c show a typical example of the evolution of (r^r 2 ), (pi,p2) 
and (V, a), respectively. Units, s.t. U — 1 and a — 1, have been chosen. We have 
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set the initial conditions for (r^r 2 ), p\, a and V. The initial condition for p2 
has been decided using the constraint equation eq.(21). We can observe the same 
asymptotic behavior of the system which arises irrespective of the initial conditions, 
due to the back-reaction: The back-reaction drives the system into the direction 
corresponding to a thinner torus, i.e. r 2 — > while r 1 — > finite. At the same time, 
the 2-volume V asymptotically approaches zero. We find out that this behavior is 
universal by setting various generic initial conditions. This universal behavior can 
also be understood by investigating the qualitative characteristics of eqs.(21)-(24), 
which shall be done in the next sub-section. 

We should also note a special class of trajectories characterized by the initial 
condition r 1 = n/2 (n : integer), p\ = 0* The (r 1 , r 2 )-trajectory becomes parallel 
to the r 2 -axis and (pi,p2)-trajectory is on the j>2-axis. Depending on whether 
P2 > or p2 > 0, r 2 tends to oo or 0, respectively. In any case, the shape of the 
torus becomes thinner and thinner as it evolves. (Note the modular invariance of 
the system.) 

(c) The asymptotic behavior of the system 

We can understand the universal behavior of the system by looking at eqs.(21)- 
(24), and investigating the asymptotic behavior of the system as t — > oo. Key 
behaviors are 

(i) V — > 0, cr — > oo, aV n -> oo (n = 1, 2, 3, • • •)■ 

o 

(ii) {t 2 ) 2 (p± + p 2 ,) increases, at least as fast as a 2 V 2 . 
(Hi) t 2 I 0, or t 2 — > oo. 

* Because of the modular invariance of the system, the cases of r 1 =integer are equivalent 
to the case of r 1 = 0, and those of r 1 =half integer are equivalent to the case of t 1 = 1/2. 
The trajectories of the former cases are stable against perturbations, while the trajectories 
of the latter cases are unstable. This can be seen from Figures 1 — a, b along with eq.(17). 

§ Here, l y(t) increases at least as fast as x(ty or l y(t) increases at least as #(t)' means that, 
\x(t)/y(t)\ — > c, < c < oo when t — > oo. In other words, l/y(t) = 0(l/x(t)) when t — > oo. 
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(iv) pii 1 + P2T 2 increases at least as a 2 V, and -^^((f 1 ) 2 + (f 2 ) 2 ) increases at 
least as a 2 . 

Now, let us derive the above results. First of all, eq.(22 — b) can be written 
with the help of eq.(21) as 

a = aa 2 + -(r 2 fl 2 f(r)V^I 2 . (22 - b>) 

Thus, cr > 0, so that a always increases and becomes positive at some stage. 
Then, V decreases because of eq.(22 — a). Furthermore, it is easily shown that 
(a 2 V 2 ) = h(r 2 ) 3 / 2 f (t)ctV 1 / 2 > 0, so that the combination o 2 V 2 always increases 
in time. (Therefore, a 2 V increases more strongly.) Afterwards, it is easy to get 
(i), by induction. Next, from eq.(21), (t 2 ) 2 (p\+ p 2 ) increases, at least in the same 
manner as a 2 V 2 . Thus, we get (ii). Now, from eq.(24 — b) and eq.(21), 

P2 = '^Ap 2 + p 2 ) - ^a 2 V + h(T 2 f 2 ^V-V 2 , (24 - b>) 

so that p2 < (note that < 0). Furthermore, |_p"2 1 > j^v 2 V, so that 

P2 decreases faster than — g ^r-- (Note that, if r 2 ^finite, this implies a strong 
deceleration of P2-) This fact excludes the behavior r 2 — ^-finite, for, if so, p2 oc 

I ''2. \2 

^ (^2)2 (eq.(24 — a)) should tend to zero, which contradicts with the deceleration 
of p2- Therefore, r 2 always behaves as r 2 j or — > oo. Thus, we get (Hi). Next, 

/ 2 \ 2 

we see that p\r l + P2T 2 (= v (p 2 (eqs.(23 — a), (24 — a)) increases at 

least as a 2 V, with the help of eq.(21). Finally, ^((f 1 ) 2 + (f 2 ) 2 ) (oc p x r x +p 2 r 2 ) 
increases at least as a 2 V, so that -^^((t 1 ) 2 + (f 2 ) 2 ) increases at least as a 2 . Thus, 
we get (iv). 

The generic trajectories are the ones for which r 1 ^finite and r 2 — > 0, like 
Figures 2 — a,b, c. We can understand this behavior as follows: Suppose that If 1 ] 

9 • 2 

is at most comparable with |f |. Then, from (iv), we can make an estimation as ^ 
~ —a, so that r 2 rapidly approaches to (faster than exp — at since a is increasing). 
Noting that V^ 1 increases much slower than a ((i)), the combinations of the form 
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(j 2 ) n V~ m in (23 — a, b) become strong suppression factors. This is compatible with 
the assumption that \ f 1 \ is not so large compared with |f 2 |. Therefore, in (24 — 6), 
only the term proportional to p\ on the R.H.S. dominates and determines the gross 
properties of the equation, which gives rise to the universal behavior. 

There is a special class of trajectories determined by the initial condition r 1 = 
(or in general, r 1 = n/2 (integer)) and p\ = 0* Due to the property of /(r^r 2 ) 
(eq.(18)) with eqs.(23 — a,b), this implies that r 1 = (or = n/2), p\ = 0, i.e. 
the trajectory of (r^r 2 ) and (pi,P2) form a line-segment on (or parallel to ) the 
r 2 -axis and ^2-axis, respectively. Combining (iv) with p\ = 0, we see that pit 1 
always increases. It means that any (r 1 , r 2 ) -trajectory which is parallel to the 
r 2 -axis has no turning point, and that r 2 tends to or oo, depending on the initial 
condition. Furthermore, combining again (iv) with p\ = and r^constant, we 

' 2 o 

see that ±<r, so that r 2 approaches rapidly to oo or (faster than exp±at 
since a is increasing). 

As is noted previously, our treatment is based on the adiabatic approximation. 
Thus, the results should always be taken with a caveat. In general, when instability 
is observed in the adiabatic treatment, it implies the unstable tendency of the 
system and it suggests the necessity of a further investigation beyond the adiabatic 
approximation, rather than just neglecting the resultant instability. Furthermore, 
in the present case, there are good reasons to regard the unstable behavior as a real 
one. First, as investigated above, the universal asymptotic behavior of the generic 
trajectories implies that f 1 — > 0, r 2 — > ( and V ~ o(er)) although pi — > — oo. This 
is because f 1 oc ^" y Pl , f 2 oc ^ T y P2 , and r 2 becomes a strong suppression (stronger 
than expert), while 1/V is at most ~ a. Thus, the adiabatic treatment for r 1 and 
t 2 becomes better and better as r 2 0: cj^/cj^ 2 ~ A^/A^ 3 / 2 ~ (r 2 ) 3 / 2 • -^y = 
(r 2 ) 1 / 2 • (^) — > (see eq.(5)). Furthermore, K does not harm the adiabatic 
treatment because of the conformal invariance of the matter field, as has already 
been discussed previously (§§2 — (b), after eq.(8)). See Figures 3 — a, b, c. (On the 

* The remarks for the last paragraph of §§3 — b apply here, too. See the footnote there. 
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other hand, we should also note that the special class of trajectories characterized 
by r 1 = n/2 (n: integer), r 2 — > oo, is not appropriate for the adiabatic treatment: 
By (iv), r 2 tends to infinity even stronger than expat. However, because of the 
modular invariance, the trajectories for which r 1 =constant and r 2 | give the 
good information of the class of these trajectories.) 

Another support for the present result comes from the consideration of the case 
of the negative cosmological constant without matter field. It is straightforward to 
introduce the A-term [24] (see eq.(16)): 

aH = ^^(p 2 +p 2 )-^a 2 a 2 V-aAV = , (21A) 
V = -aaV , (22A - a) 

- = V + ^(p?+p|)+A , (22A-6) 

f^^r 2 ) 2 ^ , (23A-a) 

p\ = , (23A - b) 

r 2 = ^(r 2 ) 2 P2 , (24A-a) 

P2 = -yr 2 (p 2 +p 2 ). (24A-6) 

Here, —A corresponds to the cosmological constant (A > 0). Because of the nega- 
tivity of the last term in eq.(2lA), the same kind of evolution for (r 1 , r 2 ) as in the 
case of the matter field is observed. (It is also notable that (21A)-(24A — b) can 
be solved analytically [24].) It strongly suggests that the instability is indepen- 
dent of the adiabatic treatment. At the same time, we should note the essential 
difference between our case and the case of the negative cosmological constant. Es- 
pecially, the difference between (23 — b) and (23 A — b) is prominent. Furthermore, 
/i(t 2 ) 3 / 2 /(t)\/~ 3 / 2 (y^h corresponds to A comparing (21) with (21A)) depends 
on (r^r 2 ) and V, which causes a highly non-trivial evolution. 
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4. The effective action for the modular degrees of freedom 



(a) Partition function 

We have treated so far the back-reaction of the quantum field on the modular 
degrees of freedom, in the sense that the semiclassical Einstein equation, eq.(l), 
has been solved, with < T a p > on the right-hand side being calculated in the 
background spacetime. We can handle the same problem in a more systematic 
manner by the path-integral approach. The significance of this investigation is as 
follows: 

First, we know that we can derive eq.(l) formally, by taking the first variation 
of the phase w.r.t. g a p in the in-in path-integral expression for g a p and ip [5] [6] [7]. 
However, when we discuss the semiclassical gravity in more detail, it is preferable 
to take into account the effects coming from the path-integral measure of g a p. 
Since we cannot fix the measure in a reasonable manner, we usually do not discuss 
much about this effect. Fortunately, our model is simple enough to investigate the 
measure to a great extent, by making use of the techniques developed in string 
theories [15] [16]. 

Second, regarding this problem, we expect that the measure in the original 
phase space, J[dh a i)dTT ab dNdN a ], should reduce to the standard canonical mea- 
sure in terms of the reduced phase space variables, J[dr A dpA][dVda][dN], after 
gauge-fixing, according to the general theory of the path integral for the lst-class 
constrained systems [25] . Analyzing this reduction process in detail for the case of 
our model is highly non-trivial and helpful for deeper understanding of the path- 
integral approach for quantum gravity [26]. 

Third, furthermore, our model also becomes a test candidate for another fun- 
damental problem: The validity of the minisuperspace approach in quantum cos- 
mology. It is essential in our reduction procedure that the condition of iV =const 
on E is compatible with the equations of motion (§§3 — b). In the context of 
quantum cosmology, it can correspond to the minisuperspace approach: We often 
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impose the special form on metrics, which is compatible with the equations of mo- 
tion, and quantize them within this sub-class of metrics, for tractability. Then, 
a fundamental question arises as to whether this approximate treatment reflects 
faithfully the main features of the full-quantized system. The results may depend 
on which space is chosen as the starting whole phase space, viz. whether we start 
from the full phase space (full quantization) or from its sub-space (minisuperspace 
quantization). In the former case, it is expected that some extra factor emerges in 
the measure, since in this case, the condition N =const on E should be treated as 
an extra constraint, rather than just an ansatz. If so, this extra factor can have 
some influence on the semiclassical evolution of the system. The similar effect can 
arise from our assumption of the spatial homogeneity of our torus model (§§2 — a). 
Our model is suitable for the detailed analysis of this fundamental problem. In the 
present paper, however, we restrict ourselves to the treatment a la minisuperspace 
models, which itself is one consistent treatment. 

Fourth, when we need to investigate validity conditions of the semiclassical 
treatment described by eq.(l), then, we have to study the second variation of the 
effective action VK[V^, t + ; VL, t_] [7]. Thus, we need to estimate W[V+, r + ; V_, r_] 
using the in-in path-integral formalism. 

We first discuss within the framework of the standard in-out path-integral 
formalism [8] and later generalize it to the in-in formalism. In this subsection, 
we shall derive the expression for the partition function Z in terms of the reduced 
phase space variables. In the next subsection, we shall estimate the effective action 
for matter, W[V+, t+; V_, t_]. 

The partition function in our case is given by 



where, in the last line, we understand that the matter degrees ip have been inte- 
grated out and suitable vacuum expectation values have appeared in TC and TC a 



Z = M I [dh ab dn ab dNdN a ] [dip] exp % \ {7T ab h ab + p^tp -NH- N a H a ) 



^ a hah — NH — N a H a ) 




(25) 



23 



(e.g. T a pn a n' — >< T a p > n a n ' ). (See the next subsection for more explicit 
discussions.) 

Integrating over the multiplier N a is equivalent to inserting S(TC a ) and setting 
N a to be an arbitrary value if needed* Let us set N a = 0: 

Z = Af J [dh ab d7i ab dN]5(H a ) exp i J (n ab h ab - NH) . (26) 

The action is invariant under the time-reparametrization and Diff(E) (the diffeo- 
morphism on the spatial surface £). Now, the gauge-fixing is needed to make this 
expression meaningful. The gauge-fixing condition for Diff(E) which is directly 
connected to our classical treatment in §3 is, 

h ab -Vh ab = , (27) 

where h ao is given in (2 — 6). 

At this stage, we need to fix our general attitude for the treatment of our model. 
Any 2-dimensional metric h ab is conformally flat [15] [16], and the conformal factor 
V is a function of spatial coordinates (as well as a time parameter t) in general, V = 
Vft,^ 1 ,^ 2 ). Here, furthermore, we set a further restriction to construct a tractable 
model, which we have investigated in the previous sections: we restrict the class of 
spatial metrics h ab to the one in which V becomes spatially constant, V = V(t). 
At the same time, the lapse function is restricted to N — N(t). Both of these 
ansatz are compatible with the classical equations of motion. Such restrictions on 
the class of the path-integral variables correspond to the minisuperspace models 
in quantum cosmology. (See §5 for more discussions on this point.) 

The treatment for the time-reparametrization invariance is well-investigated 
[27]. The final result is neat: Introducing the physical time T = f dt N(t), one 

* This situation is parallel to the case of QED. In the latter case, the term AodivE appears 
in the action. One can set A = if needed, provided that S(divE) is inserted in the 
integrand. 
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computes a transition amplitude from time to time T. Then, integrate over 
the result with respect to T [27]. Here, we shall not do it explicitly, since we are 
mainly interested in the semiclassical evolution of the system. We understand that 
we follow the above procedure whenever needed. 

Then, 

Z=M J [dVdv' a d 2 r}[dh ab d7r ab dN}5(h ab - Vh ab )A F p5((Pl%) a ) expiS 

= Mj [dVdv' a d 2 r][d^ b da]J[dN]A Fp lhab=vh J((P!*) a ) ewiS ]hab=vhab , 

where S = J ir ab h ab — NH. Note that h ab = V(t)h ab corresponds to choosing the 
York's time-slicing, K = TT a a /V =const w.r.t. the spatial coordinates [21] (see 
§§3 - b). Thus, only the traceless part of vr a& , Tc ab = T\ ab - \nh ab = K ab V, has 
remained in argument of the delta-function in the last line above. Accordingly, the 
change of the integral variables has been performed and J is the 

Jacobian factor associated with this change. Employing the method in Appendix 
-B, J can be determined as follows: A natural diffeo-invariant inner-product* for 
5n ab is (Sir, 5ir) = J d 2 ^Vhh ac h bd 5n ab 57T cd . Substituting n ab = fc ab + \h ab aV, we 
get (5tc,5tt) = (5n,5n) + ±V 2 (5a,5a), where (5a, 5a) = J d 2 ^Vh(5a) 2 . Thus, 
1 = J J dTT ab daexp—(5ir,5-K), so that J = V up to an unimportant numerical 
factor. 

The Faddeev-Popov determinant App in our case is equivalent to the Jacobian 
associated with the change of the integral variables from h ab to (V, v' a , (V 1 ,! -2 )), 
where v' a KerP\. Thus, we can employ the method in Appendix B again to 
determine App: From eq. (Al), 

1 1 2 = \ \5<ph ab + (PW)ab + T Aab 5r A \\ 2 

= 4(50, 50) + (v', PlPW) + (T A , T B )5r A 5r B . 



* An appropriate power of a := Zpi an ck should be multiplied to formulae in order to adjust 
physical dimensions like eq.(A3). It is easy and not significant for the present discussions, 
so we omit the factor. 
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Then, 



1 = A FP J d(f)dv'd 2 T&cp-\\Sh\\ 2 = AFpidet'PlP^-^det-^^TB) 
Thus, 

A FP = (det'PlP^det^iTA, T B ) ■ 

Thus, 



Z = Vol Difh M J [dVd 2 r][dn ab da][dN] 

(i t \ 1//2 
* (i/2 ( T ^>™« p i , *>>*p' s i/...^„, . 

(28) 

where, {x a }a=i,2 is the basis for KerPi, a space of conformal Killing vectors* 

Let us investigate the factor <5((P^7f) a ). According to eq. (CI) in Appendix C 
(A = P\, x = n ab , f(x) = expiS and {^ A }a=i,2 are the zero-modes for P-J), 



/ 



[dn ab ] S((Pln) a ) exp iS[h ab = Vh ab , ir ab , a, N] 

r dpt l / 2 (<n A r ( 29 ) 

= / [ d 2 p ] ^ I exp i / dt( PA r A + aV- NH) . 

J det'P{ J 

Here, in the last line, the non-zero-mode components of Tx ab have been set to be zero 
according to the formula (CI). This is equivalent to substituting ix ab = K ab V = 
J2aPA^ ab V into the action. Therefore, this is the path-integral version of the 
procedure of solving the momentum constraint in §§3 — b. 



* Any element in Diffo (diffeomorphism on X homotopic to 1) is associated with a vec- 
tor v a , which can be decomposed as v a = v' a + A Q x": where v' a KerP\. Noting the 
argument in Appendix A, J[dv a ] = J[dv' a ] d 2 \ det 1 / 2 (x a , X^)> which means Volr>ifj a = 
(f[dv' a }) ■ VoIkctPx- Thus, by factorizing (j[dv' a }) = VolDiff /VolKerP! from the path- 
integral, the factor det -1 ' 2 (x"S X^) appears. Here, the factor (J cPA) -1 is absorbed into 
the normalization Af. 
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Thus, 



J det' llz PiPl det 1 / 2 (x a ,X fj ) 

xexp? J dt(p A r A + aV - NH) . 

(30) 

Here, ciet'P/ = (det' P^l) 1 / 2 has been used* 

Now, we choose {Ta}a=i,2 an d {^ A }a=i,2 as (T4, = S B (see §§3 — b), 
so that det 1 / 2 (T A ,T B )det 1 / 2 (ty A ,ty B ) = 1. For our case of a locally flat torus, 
{x a }a=l,2 can be chosen as x\ = (1)0) an d Xa = (0, 1), without inducing any 
critical point as vector fields. Then, <iet 1 / 2 (x Q , X 13 ) — 1- 

1 /2 f 1 /2 i" 

Finally, det' P\P\ and det' P\P\ should be estimated. The map P\ is a 

map from a space of 2- vector fields to a space of 2nd rank, symmetric and traceless 

tensor fields, and the map p\ is a map from the latter space to the former space. 

Note that each of the spaces can be represented as a 2-component vector fields. 

Now, it is convenient to use the complex coordinates (z,z), with respect to which 

both Pi and p\ become diagonal [15]. Let z = x + iy,z = x — iy. Then, the line 

element becomes (e^ := V), ds 2 = e^h^d^d^ 2 = e^{dx 2 + dy 2 ) = e^dz dz, so 
( Je*\ 

that hat = I 1 1 . (The suffix (z, z) is for the explicit indication of the 

o 

coordinates employed.) The following arguments are valid for a general spatial 



* This equality can be shown by estimating an integral I = J dw lab exp — {pIw',pIw') in 
two different manners (here, w' ab is symmetric, traceless and ^ KerP\): One way is I = 
J dw' exp -(to', Pi p/to') = {det'PxPl)- 1 / 2 , and the other way is 7 = / d(P}w')(det' P})- 1 
exp — ) = {det' P\) 1 . This change of the integral variables in the latter estima- 

tion is valid since the space of the original variables (w' ab ) is isomorphic as a vector space 
to the space of the new variables ((P/u/) a ) by the map P\. See below. 

§ We shall use the following facts: d :— d z — \(p x — id y ) and d := d z — \{d x + id y ); 
v = (v 1 , v 2 )( x ^ = (v 1 + iv 2 , v 1 — iw 2 )( Zj 2), i.e. v z = v 1 + iv 2 , v z = v 1 — iv 2 = v z (y 1 , v 2 e R); 
Let T ab be symmetric and traceless, and let its components in (x, y)-coordinates, T 11 etc., 
are real, then {T ab ) [z ,- z) = diag (2(T n + iT 12 ), 2(T n - iT 12 )), i.e. T zz = 2(T n + iT 12 ), 
T zz — 2(T n — iT 12 ) — T zz and the other components vanish; The Christoffel symbols 
become T z z = d<j>, r§ f = 3(f> = T z z and the others vanish. 
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metric h a j, on a torus, so that we shall discuss in general terms. Only at the final 
stage (eq.(34) below), we set the condition that = InV = spatially constant. 

Now, both P\ and p| can be regarded as a map from a 2-component field to 
another 2-component filed: 

Pi : t (v z ,v z ) i — > ^{Pivy^iPivf"), P\ : \w zz ,w™) i — ► ' ((Pjiu)*, (P^iu)*), 
where w ab is a symmetric, traceless tensor field and *(•, •) indicates the transposi- 
tion. In this sense, Pi and P| are represented as 

Pl= { 4e"*0 J ' Pl = V -2e-^de 2 * ) ' (31) 

Thus, P|Pi : t (v z ,v z ) i — ► * ({pIPiv) z , (p/Pi)*) is represented as 

PlPl= l^ -8e-^9j = V 2A + ( 2 )pJ ' (32) 

where A = —D a D a and ( 2 -*P are the Laplacian and the scalar curvature, respec- 
tively, defined by the covariant derivative (D a ) w.r. t. e^Kb. Similarly, PiP\ : 
t (w zz , w zz ) i — > 1 ( i (PiPlw) zz , (PiP\w) zz ^j is represented as 



PiP\ 



' Se-tde-Wd \ /2A-2(% 

-8e-^9e- 2< ^a / ~ V 2A - 2 ( 2 )p 



(33) 
Therefore, 

det' l / 2 p\Pi = det'{2A + ^R) , det' 1/2 PiP\ = det'{2A - 2^R) . 
In our model of locally flat tori (0 = \nV = spatially constant), thus, 
det' 1/2 P\Pi = dei! xl2 PiP\ = det'(2A) . 
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Finally, we obtain 

Z = M J [dr A dp A ] [dVda] [dN] exp i J ( PA f A + aV - NH) . (34) 

The integral region for r^r 2 should be understood as on the moduli space, 
M. g =\: As is indicated in eq.(28), Diffo(E) (the diffeomorphism group on E no- 
motopic to 1) has been factorized from the path-integral. What is really needed to 
be factorized is the whole diffeomorphism group on E, Dif '/(E). Note that [15] [16] 

M g ~ Riem(Y,)/Weyl x Diff{T) ~ (iftem(E)/WeyZ x Dif f Q {T)) /MCG 
~ H + /PSL{2, Z) ~ D(H+)/ ~ . 

Here, MCG := Dif f(Y)/ Dif f Q (Y) is the mapping-class group for E, and MCG ~ 
PSL(2, Z) for E ~ T 2 (i.e. a group of 2 x 2 unimodular matrix with integer 
elements, modulo sign). D(H+) is the fundamental region in H+ (upper half-plane) 
w.r.t. the action of PSL(2,7i) (e.g. the Dirichlet region D = {z E H+\\Rez\ < 
1/2, \z\ > 1}) and / ~ indicates the identification (t^t 2 ) ~ — (r^r 2 ) on the 
boundary of D [15] [16]. Thus, the integral region for (r 1 ,^) in eq.(34) should be 
understood as over A4 g =i rather than over H+, considering that we have factorized 
the volume of the mapping-class group MCG ~ PSL(2,Z) as well as Z)i//o(E) 
from the path-integral. 

If we integrate out the momenta pa and a in eq.(34), we get 
^ = ^/[^][^][^^</*^W{^^((0 2 + (r 1 ) 2 ) (35) 

1 v 2 + h{T 2 fl 2 f{T)v- 1 ' 2 } . 
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Note that the kinetic term for (r 1 ^ 2 ) in the action is proportional to Gabt A t b 
and the last term in the action is proportional to the non-holomorphic Eisenstein 
series G(r, 3/2) (see below eq.(18)). Thus, Z is modular invariant since both the 
measure and the action are modular invariant. 
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(b) Estimation of the functional determinant for the matter 

Now, we estimate the path- integral for the matter ip in eq.(25). Our aim is to 
obtain the effective action of the form W[<f>, r 1 (-), r 2 (-)] by integrating out quantum 
fluctuations of the matter. Generalizing the framework to the in-in formalism and 
getting r + ; </>_, r_], one can discuss the validity conditions for the semiclas- 

sical treatment [7], eq.(l). At this stage, the peculiarity of the system including 
gravity is prominent. In the standard treatment of a dissipative system, like a 
quantum Brownian motion [17], the interaction between the sub-system and the 
environment is described by a weak, linear coupling. In our case, however, there 
is no such interaction term between gravity (analogous to the sub-system) and 
matter (analogous to the environment). Rather, the interaction is bilinear in tp 
and non-linear in (r^r 2 ) and 0, as is seen from eq.(7). Thus, it requires a new 
treatment for a deeper analysis. Here, we should be content with only a rough es- 
timation of the effect of the nonlinear coupling. We want to estimate the partition 
function for the matter, 

Z^ = J [#] exp ~ J V(-<9 2 + l -R)^ = Det ~ 1/2 i^-(-d 2 + \r)} 
= exp~W[r(-)] • 

Here, " ~ " denotes the Riemannian signature quantity. We calculate using the met- 
ric 9a/3 — (1> Vhab) with eq.(2 — b). It is difficult to estimate the above functional 
determinant exactly for a general function (r 1 (-), r 2 (-)) and V(-). From the view- 
point of the quantum dissipative system, this difficulty comes from the peculiarity 
of the interaction between gravity and matter. As discussed in the beginning of 
§§2 — b, we treat the back-reaction problem in the sense that we investigate the 
modification of the background geometry due to matter, i.e. due to < T a p > 
calculated on the background spacetime. We have chosen as a background, a flat 
spacetime. Thus, for the lowest order approximation, we treat r 1 , r 2 and V as con- 
stants, so that we can set R — 0. This treatment corresponds to the lowest order 
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estimation of the functional form of the effective potential in standard quantum 
field theory [28]. 

Thus, we need to estimate the determinant of the operator 

.._ a-v B ^_^v (g2g + v _ lhabaA) _ 



2nh 2nh 

where % and a 2 have been inserted for the convenience of recovering a formula for 
pseudo-Riemannian signature. Now, we need to solve the heat equation 



A P=~lsP » 

lim slo p(x, y, s) = <5 (3) (x - y) . 

Here, x := (x° = t,£\£ 2 ). Taking care of the periodicity in space, the solution is 
given by 

h ^ 



£ ex p -^v~s {{x ° ~ x °' )2 + vhab{i ~ i>+ nm ~^' + n)bh 



X 

ni,n 2 



especially, 

3/2 



x 7 ni,n 2 

where (n,n) := h a bn a n b = -pi(n\ + 2r 1 nin2 + [rpn 2 ,). Thus, the ^-function associ- 
ated with A is [27], 

00 

Ca(z) = f^J dss z - l Trp(s) 





oo 



£ / ^ { 

n, no J V 



r(z) ^ 7 V 2a 2 Fs / F 2a 2 



nm 2 g 



ex P-7T7^>> n ) 
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where Q = J d^Xy/g and a transformation of variable s (x := ^-(n, n)s" 1 ) has 
been done to get the formula in the last line from the middle line. Noting that 
£\ z=0 ( j^i) CZ ~ 1 )) = ^T C ~ l for VC7 when C is independent of z, we get 



mri2 



Thus, we get 



W = hnDetA = ~£ A (0) 



TlVt , 2\3/2 ^ tV7 \ 



4vra 3 \/ 3 / 2 ^ (nf + 2r 1 mn 2 + |r| 2 n|) 3 / 2 



To recover W for the pseudo-Riemannian signature, we replace H — > ift, a — > ia (no 
change in f2, dx°dx 1 dx 2 dx 3 <-> dx°dx 1 dx 2 dx 3 ). This replacement comes from the 
comparison between = i^lnDet (ffl( _ ^ 2 )) and = |lnDet (^[(-<9 2 )j. 
Thus, 

WV r 2 l = — fr 2 ) 3 / 2 V - (37 - 6) 

1 ' J 4^1/3/2^^ Z^ (n 2 + 2r l nin2 + | r |2 n 2 )3 /2 " ^' ^ 

Since we have used the expectation value of the energy-momentum tensor for 
the matter, < T a p >, to couple with gravity (eq.(13 — b) or eq.(16)), we need to 
use the in-in path-integral formalism, rather than the standard in-out formalism 
[5] [6] [7]. Then the matter part of the action (pseudo-Riemannian) (see eq.(7)) 
should be reinterpreted as, 

c 

= -\J (9 aP d a ^ + \m 2 )^g + \j (g a0 d a ^ + l -R^ 2 )^g , 
+ 

where " c " stands for the closed-time contour and " + " and " — " stand for, 
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respectively, the +-branch and the — branch of the time-contour. Then, 




+ 




1>- 



Since + and — components are separated completely, it is enough to look at only 
the +- sector (or — sector). 

Now, let us investigate the effective action, Spr 1 , r 2 , V; N] = S^r 1 , t 2 , V; N] 
+W[t 1 , t 2 , V; N], where S g [(j), r 1 , r 2 ; N] is the reduced action for gravity in terms 
of the configuration variables and W^r 1 , r 2 , V; N] is given by eq.(37 — b). The 
effective action S^r 1 , r 2 , V; N] is what has appeared in the exponent in eq.(35). It 
should be noted that the first variations of S[(f>, t^t 2 ] w.r.t. N, V and (r^r 2 ) 
reproduce exactly eqs.(21 — 24). This result shows the following two points. 

First, our approximation for the estimation of Det~ x l 2 (^(— d 2 + |-R)j, treat- 
ing r^r 2 and V as if they were constants so that R = 0, corresponds to the 
approximation used to solve the semiclassical Einstein equation, eq.(l). Namely, 
< T a p >, calculated on a flat background, is used in eq.(l) to estimate the devi- 
ation from the original background geometry. As is discussed at the beginning of 
§§2 — b, the latter approximation has been implemented for the tractability of the 
problem, at the expense of the self-consistency of eq.(l). Such an approximation 
is what is usually meant by the term "back-reaction" , and this may be the best we 
can do in practice. 

Second, regarding the path integral expressions in the Lagrangian formalism, 
like eq.(35): We can reproduce eq.(l) (or equivalently, eqs. (22-24)) from the phase 
part S g + W in the partition function Z with the matter part integrated, and 
without taking care of the contributions from the measure for V and (t^t 2 ) (see 
eq.(35)). However, we now know explicitly the non-trivial path-integral measure 
for V and (r 1 , r 2 ) as is shown in eq.(35). There should be 0(h) correction to eq.(l) 
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coming from the path-integral measure for g a p and this correction will cause a non- 
trivial correction to the dynamics of g a p. We shall come back to this point in the 
next section. 

5. Discussions 

In this paper, we have investigated the semiclassical dynamics of the topological 
degrees of freedom, (r^r 2 ), which has been seldom discussed so far. By reducing 
the spacetime dimension to 3, we could concentrate on the study of a finite number 
of topological modes and we could describe the back-reaction effect from matter 
to topological modes, explicitly. We observed a non-trivial dynamics caused by 
the back-reaction. The back-reaction makes the toroidal universe unstable: The 
shape of the torus becomes thinner and thinner, while its total 2-volume becomes 
smaller and smaller. These are universal behaviors of the system independent of 
the initial conditions, which is justified by the asymptotic analysis of the set of 
dynamical equations. This observation implies the importance of the investigation 
of topological aspects for a deeper understanding of quantum gravity. Moreover, 
we could fix the path-integral measure for (r^r 2 ) and V and observe that the 
partition function is expressed in terms of the canonical variables for the reduced 
phase space with the standard Liouville measure. 

Let us note a few points regarding the path-integral measure. 

We obtained the path-integral expression on the reduced phase space with 
the Liouville measure (eq.(34)), while the path integral on the configuration space 
requires a non-trivial measure (eq.(35)). Indeed, the combination ^Jw is essential 
to make the partition function modular invariant. In our model, the semiclassical 
Einstein equation, eq.(l), corresponds to eqs.(21) — (24), and they are derived 
from the variation of the exponent in eq.(34) or eq.(35). It means that, from 
the viewpoint of the Lagrangian formalism, the semiclassical Einstein equation is 
derived from the variation of the phase part in the partition function, with the 
measure factor untouched. Thus, the measure factor gives the 0(ft)-correction 
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to eq.(l). In our model, the term J dtN(t)h(2lnr 2 + 1/2 In V) can be added to 
the action as a correction. (Note the time reparametrization invariance implied 
in eq.(34).) Then, it is a non-trivial question worth while to investigate which 
is better as the semiclassical description, the semiclassical Einstein equation in 
terms of the canonical variables (eqs.(21)-(24) in our case), or the same in terms of 
the configuration variables with suitable corrections originating from the measure. 
If we perform the path integral exactly, both the canonical and the Lagrangian 
formalisms will give equivalent results, but they will not be equivalent within the 
accuracy of the semiclassical approximation. 

Another important problem is linked with the validity of the minisuperspace 
treatment. We have investigated the homogeneous model, which is equivalent to 
assuming N = N(t), V = V(t) (see the discussion in §§4 — a, below eq.(27)). 
We can set this ansatz since it is compatible with the dynamics. This treatment 
corresponds to the minisuperspace approach in quantum cosmology. Though such 
a treatment is completely self-consistent, it is important to question to what ex- 
tent such a treatment reflects the original full quantum theory faithfully. From 
the viewpoint of the original full system, the restrictions are regarded as extra 
constraints on the phase space. These constraints can modify the path integral 
measure for the reduced variables (minisuperspace variables). Since this problem 
is a fundamental one, it should be investigated separately. Our model may be a 
good test candidate to investigate this point in detail. 
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APPENDIX 



A. Brief summary on the moduli space 

We give here a concise summary on the moduli space just for fixing the ter- 
minology and notations used in §3 and §4. See e.g. [15], [16] for more detailed 
information. 

Let S be a 2-dimensional, compact, closed, orientable manifold with genus g. 
The moduli space M g of £ is defined as M g ~ Riem(T) / W eyl x Dif 7(E), where 
Riem(Tj) is a space of all Riemannian metrics on S, Weyl is for the Weyl group and 
Dif f (11) is for the diffeomorphism group on £. The universal covering space of 
Ai g is called the Teichmuller space. Now, the tangent space of M.g, T(Ai g ), can be 
investigated as follows: Any variation of the spatial metric 5h ab e T(Riem(T>)), can 
be decomposed into the trace part and the traceless part, the latter is furthermore 
decomposed into the diffeomorphism 5r>h ab , and the moduli deformation SMh ab ; 

dh a b = $whab + &Dh a b + &Mh a b 

= 5 w hab + (Piv)ab + T Aab 5r A , (Al) 

where 

SwKb = HKb for 350 , (A2 - a) 



(P 1 v) ab = D a v b + D b v a -D c v c h ab for 3v a , (A2 - b) 

dh ab _ 1 d dhcd 
dr A 2 dr A 



T Aab = ~77~~T - ~ h ^ZA h ab ■ (A2 - c) 



Here {t a }a=i .. .,dimnM g are the Teichmuller parameters specifying a point in M.g- 
A natural inner-product on T(Riem(Y)) is introduced as 

(A,B):=± J d 2 xVh h ac h bd A ab B cd for VA ab , \/B ab e T(Riem(E)) , (A3) 



where a is the Planck length, inserted to adjust the physical dimension. Then, the 
tangent space of the moduli space, T(M g ), can be characterized by the set of all 
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symmetric, traceless (covariant) tensors which are perpendicular to T(Dif f(H)) 
w.r.t. the inner-product (A3), the latter condition being equivalent to the condition 

(pjw) a = -2D b w ab = , {AA) 

for w G T*(M. g ). Thus, dim^M. g = dim^iT* (M. g ) = dim-^K er P\\ which is 
known as = 0, = 2, and = 6g — 6 for g — 0, g — 1, and g > 2, respectively. It is 
also known that dim^KerP\ — dim^KerPY = 6 — 6g (Riemann-Roch theorem). 
For the case of a torus (g = 1), then, dim^M = 2 and dim^iKerPi = 2. Thus, two 
Teichmuller parameters (r^r 2 ) are needed to describe the modular deformations 
^Afhab £ T(M. g= i)i an d two independent vectors {x a }a=i,2 are needed as the basis 
of KerP\. 

Let {TAab}A=i,2,-,dim R M g be the basis of T(M g ), and {^ Aab }A=i,2,-dim R M g 
be the basis of T*(M. g ). They can be chosen to satisfy (\P , Tg) = 5 b- Then, 
they define a metric on M. g =\ (the Weil-Peterson metric), induced from the inner- 
product eq.(A3) on T(Riem(Y>)), 

Gab = (Ta,Tb) , 

(A5) 

g AB = (\$ A j \£ B ) = inverse matrix of Qab ■ 

B. The Jacobian associated with a change of integral variables. 

Let us note a convenient method to specify the Jacobian associated with a 
change of integral variables. (See e.g. [15].) 

If a line element ds is given on a space of integral variables (X A , A = 1, 2, • • • , n) 
as ds 2 = G ABdX A dX B =: (dX, dX), then d n X\fdetG is a natural integral mea- 
sure, where V det G takes care of the Jacobian factor. Suppose we change the 
variables from X A to X A , then d n XW det G' is the corresponding integral mea- 
sure for the new variables. Now, a convenient way to find out the expression for 
Vdet G' is 

1) Express 5X A in terms of 5X A \ 5X A = &5X A> . 
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2) Then express (SX, SX) in terms of SX', (SX, 5X) = §§w§$w(SX A ' , 5X B '). 
(This should be equivalent to Ga>B'$X a '5X b '.) 

3) Then determine the Jacobian J by setting 1 = J J d n 5X' exp — (5X, 5X), 
since this should be equivalent to 1 = J(det G' /n)~ l l 2 . (The factor 7r is 
usually unimportant and omitted.) 



C. A formula for the delta-function. 

Let us derive a formula which modifies an integral including 5 (Ax) into a more 
practical form. Here, A is a linear operator possibly with zero-modes. 



Let us consider an integral, I = f dx 5(Ax)f(x). Let } (A = 1, 2, 



m 



dimKerA) be the zero-modes for A. Then, any element a; of a vector space V can 
be decomposed as x = X + J2 a Pa^ A , where X G V / KerA. Now we change the 
integral variables from x to (X, pa)- Then, (x, x) = (X, X) + (^ A , ^ b )paPb, where 
(•, •) is a suitable inner-product, which is assumed to be given. Thus, according to 
Appendix B, the associated Jacobian J becomes J = det 1 / 2 (^ A , ^ B ). Thus, 

1 = J dXdpdet 1 / 2 (^ A ,^ B )5(AX)f(X, p) 
= J dpdet l l 2 (^ A ^ B )(det'A)- l f(X = 0, p) , 

where an equality, 5(AX) = (det'A)~ 1 5(X) when X £ V/KerA, has been used in 
the last line. (This equality can be shown easily by the variable change from X to 
Y = AX). Therefore, we have obtained a formula 

Jdx5(Ax)f(x) = Jdpdet 1 / 2 ^ A ,^ B )(det'A)- 1 f(X = 6,p) . (CI) 
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Figure Captions 

Figurel-a: The plot of the function /(r 1 , r 2 ) for the range r 1 : - 1 and r 2 : 0.5 - 0.8. 
The infinite summation has been truncated at —200 and 200. 

Figurel-b: The contour plot of /(r 1 , r 2 ), with the same range and the truncation points 
as in Figure 1 — a. The lines indicate the values (from bottom to top) 30, 
28, 26, 24, 22, 20, 15 and 5. 

Figure2-a: The trajectory of (r^r 2 ) determined by eqs.(21)-(24). The infinite summa- 
tion in the definition of /(r) has been truncated at —200 and 200. % and 
a have been set unity. The initial conditions are r 1 = 0.500, r 2 = 0.500, 
pi = 1.000, p 2 = 1.620, V = 1.000 and a = 0.000. Points A - F and Z indi- 
cate typical points on the trajectory. A: the initial point, B: a point near the 
turning point {j>2 = 0), C — F: the points for which t is near integer, Z: the 
end point of the calculation. A : (0.000, 0.500) at t = 0.000, B : (0.727, 0.888) 
at t = 0.675, C : (1.171,0.760) at t = 1.003, D : (1.551,0.161) at t = 2.057, 
E : (1.553,0.031) at t = 2.953, F : (1.549,4.8 x 10~ 3 ) at t = 4.043 and 
Z : (1.549, 1.0 x 10~ 3 ) at t = 4.983. 

Figure2-b: The trajectory of (pi,P2) determined by eqs.(21)-(24). The initial condi- 
tions are the same as in Figure 2 — a. A : (1.000, 1.620) at t — 0.000, 
B : (0.949, -3.6 x 10~ 3 ) at t = 0.675, C : (0.949, -0.603) at t = 1.003, 
D : (0.811,-5.301) at t = 2.057, £ : (-4.960,-29.22) at t = 2.953, 
F : (-9.310, -175.3) at t = 4.043, and Z : (-12.59, -820.4) at t = 4.983. 

Figure2-c: The trajectory of (V, cr) determined by eqs.(21)-(24). The initial condi- 
tions are the same as in Figure 2 - a. A : (1.000,0.000) at t = 0.000, 
B : (1.089,-0.235) at t = 0.675, C : (1.193,-0.320) at t = 1.003, D : 
(1.825, -0.452) at t = 2.057, £ : (2.774, -0.474) at t = 2.953, 
F : (4.540, -0.423) at t = 4.043 and Z : (6.597, -0.372) at t = 4.983. 

Figure3-&: The value of f 1 during the evolution shown in Figure 2 — a, 6, c. A: 0.500 at 
t = 0.000, B : 1.500 at t = 0.675, C : 1.100 at t = 1.003, D : 4.182 x 10~ 2 at 



t = 2.057, E : -9.657 x 1(T 3 at t = 2.953, F : -4.290 x 10~ 4 at t = 4.043 
and Z : -2.519 x 10~ 5 at t = 4.983. 

FigureZ-h: The value of f 2 during the evolution shown in Figure 2 — a, b, c. A : 0.810 
at t = 0.000, 5 : -5.678 x 10~ 3 at t = 0.675, C : -0.697 at t = 1.003, 
D : -0.273 at t = 2.057, E : -5.688 x 10~ 2 at t = 2.953, F : -8.078 x 10~ 3 
at t = 4.043 and Z : -1.641 x 10~ 3 at t = 4.983. 

Figure3-c: The value of F during the evolution shown in Figure 2 — a, 6, c. A : 0.000 
at * = 0.000, 5 : 0.257 at t = 0.675, C : 0.382 at t = 1.003, D : 0.825 at 
t = 2.057, E : 1.315 at i = 2.953, F : 1.922 at t = 4.043 and Z : 2.454 at 
t = 4.983. 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/95 1 1 046v 1 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/95 1 1 046v 1 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/95 1 1 046v 1 



This figure "fig2-2.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/95 1 1 046v 1 



